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On the Law of Gravitation in the Binary Systems. 

By Frank Loxley Griffin. 



§1. The Problem. 

That the Newtonian law of attraction operates in the solar system, was long 
ago shown to be very highly probable; the discovery of binary systems naturally 
raised the question as to whether the same law operates in them. Any reply to 
such a question is based, of course, upon both the observational data and certain 
fundamental assumptions expressed or implied. In this paper the assumptions 
hitherto used in demonstrating the operation of Newton's law are replaced by 
others ; incidentally a new proof with the old hypotheses is obtained. 

The observational data have to do with the apparent orbit, which is merely 
the projection of the actual orbit upon the plane tangent to the "celestial sphere"; 
but, combined with the assumption that the actual orbit is a plane curve, they 
permit the conclusion that the satellite describes an ellipse, subject to the action 
of a central force directed toward the primary, f 

These data being accepted, the operation of Newton's law can be established 
by applying a theorem of Darboux and Halphen J and employing the following 
assumptions : the force 

(I) is a function of the distance alone, 
(II) does not vanish at the center of force, and 
(III) admits no orbits except conic sections, whatever be the initial 
conditions. 

* Bead (in part) before the American Mathematical Society Uecember 28, 1906, and (supplemented) Sep- 
tember 6, 1907. 

t F. Tisserand : TraiU de Mecanique OSleste, Vol. I, p. 35. In the last analysis, this conclusion is to be 
regarded as a pure assumption. For, in the first place, observational error is considerable in comparison with 
the small angles measured; so that the observed positions often deviate widely from even the "most probable 
ellipse". And secondly, — and more important, — even if absolutely correct observations were possible, any 
number of them would, without the help of assumptions, be powerless to give any information as to the 
positions of the star at any instants save those of the observations, or to show that the orbit is a conic with 
the law of areas exactly fulfilled. 

% Comptes rendus de VAeadimie dez Sciences (Institut de France), Vol. LXXXIV, pp. 760-62, 936-38, 939-41 ; 
cf. also Tisserand, loc. cit., pp. 36-42. 
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Another demonstration uses the following assumptions with a theorem of 
Bertrand : * the force 

(I) is a function of the distance alone, and 
(I a) is a single-valued function, f 
(II) does not vanish at the centre of force, and 
(III) admits no orbits except closed curves whatever be the initial 
conditions, provided that the velocity be less than a certain 
limit. 
For the proof given in the present paper the hypotheses are : the force 

(I) is a function of the distance alone, 
(II) does not vanish at the centre of force, and 

(III) varies along the given ellipse according to a law which, if 
operating throughout the plane, would permit real motion in 
some region outside of any circle about the center of force. 
In the present paper, as in the classic discussions, the bodies are regarded as 
material points. [See § 15, however.] 

In comparing these three sets of hypotheses, it is hardly necessary to point 
out the great severity of the third in each of the classic sets [that of Darboux 
and Halphen being more severe than that of Bertrand, involving the latter] ; 
or to notice the mildness of the third in the new set, since to deny it is to admit 
either that the force varies according to different laws in different parts of the 
plane, or that a body started anywhere outside of a certain circle with real 
initial conditions would immediately disappear from space. Also (I a), while of 
the sort readily admitted in considering the forces of nature, is one of the least 
easily dispensed with in the new set, its absence necessitating an increased use 
of the remaining hypotheses and a considerably lengthened discussion. It should, 
however, be noted that (II) can be omitted from either of the classic sets, if to 
the observational data be added the statement, probably justifiable, that the 
primary star is not at the center of the ellipse described by the satellite. And 
further, — because of a theorem due to Bertrand, J — Darboux and Halphen, 

* Comptes rendm, Vol. LXXVII, pp. 849-53; of. Tisserand, loc. cit., pp. 43-48. 

j- This is tacitly used by Bertrand in the assumption that the orbit has only one apsidal angle. For a 
simple example of an orbit with two apsidal angles see a paper by the writer in the American Mathematical 
Monthly, Vol. XIV (1907), pp. 199-301. 

% Comptes rendut, Vol. LXXXIV, pp. 731-32. 



64 Geifpin: On the Law of Gravitation in the Binary Systems. 

knowing that at least one orbit is an ellipse, do not need that part of the obser- 
vational data which states that the force is central. 

The present discussion has three parts: (1) a proof of the operation of the 
Newtonian law; (2) a treatment of the trajectories admitted by the laws satis- 
fying hypotheses (I) and (III), — an extensive category; and (3) a consideration 
of various questions relating to sets of hypotheses, including a new proof with 
the hypotheses of Darboux and Halphen. 

The demonstration in part (1) starts from the fact that the force required 
for the description of the given ellipse as a central orbit varies in a manner 
dependent upon the position of the center of force in the plane. From the polar 
equation of the ellipse, the pole being at any point, an expression is obtainable 
for the required force as a function of the distance alone. An examination of 
the equation of the ellipse shows that if the center of force be elsewhere than on 
the major axis, hypothesis (III) is contradicted; and if the center lies elsewhere 
than at a focus, (II) is not satisfied. 

It may be noted at this point that the center of force must lie within the 
given ellipse, since a closed curve can not be described as a central orbit with a 
center of force exterior to, or upon, the curve. * 

§ 2. Lemma Concerning Imaginary Forces. 

For later reference it is well to observe that if the resultant force — central 
or not — be imaginary in any portion of space, real motion is impossible in that 
region. For, suppose a particle (x, y, z) to move there, that is, let 

* = $i{*), V = H*\ a = 4>b(0, 
where the <j^(tf) are real functions. From the definition of a derivative it follows 

that t? , J and tt are real functions of t, as are also the second derivatives, or 
at ' at at ' 

component accelerations. Henoe, at every instant, the resultant acceleration 
and force are real. This contradicts the hypothesis, so that the supposition of 
real motion must be abandoned. 



* E. J. Routh : Dynamics of a Particle, pp. 
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§ 3. Polar Equation of the Given Ellipse. 

Let the center of force be any point within the given ellipse but not on a 
principal axis. Let coordinate axes coinciding with the principal axes be so 
chosen that the center of force (x lt y x ) lies in the third quadrant. Then the 
equation of the ellipse is 

( 1 — <?V + y 2 = a\\ — <?), (1) 

where a and e are respectively the major semi-axis and the eccentricity. Now 
from any point (x x , y x ) of the third quadrant there is a normal to the ellipse 
(in the first quadrant), where the radius vector is a maximum, and where, 
if % denote the angle between the X-axis (positively directed) and the normal 

7E 

(outward directed), 0<>l< -jr-. Let p x be the distance from the foot of the 

normal (x , y Q ) to (x lt y x ), and p a be the normal length to the F-axis. Then, 
by hypothesis, p!>p 2 . 

Let the center of force be selected now as pole, and the coordinates of any 
point of the ellipse be u and 0, where u is the reciprocal of the radius vector p, 
and is the longitude measured from a line parallel to OX. Then, since 

x = x x + p cos 0, y — yi + p sin 0, 
equation (1) becomes 

(1 — e 2 ) (cos d + Xj, u) 2 + (sin + y x uf — a 2 (l — <?)u % . (2) 

The elimination of sin gives 

(a cos 2 + (3 cos + y) 2 = 6 2 (1 — cos 8 0), (3) 

where 

a = — e 2 , pz=2x 1 u(l — e 2 ), y = 1 + tf x u 2 + {x\ — a 2 ) (1 — e> 2 , 

and 

h = 2y x u. 

Re-arrangement of (3) yields the desired form : 

cos 4 -f 45 cos 3 + 60 cos 2 + 4Z> cos + E = 0, (4) 



where 

/? """ 2a ' v "" 6a 2 ' " ~ 2a s ' " ~ a: 
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§4. $ as a Function of u. 
It is necessary to examine the solutions of the quartic equation (4).* Now 

this equation is satisfied by cos 6 = cos "K if u = w x = — . Since u is a mini- 
mum at 6 = \ there are not only values of greater than X but also values of 6 
less than X, for which u has any value slightly greater than u x . Hence equation 
(4) defines cos as a multiple-valued function of u ; that is, cos 6 is one function 
of u one side of the maximum normal and another function the other side. 
But cos 6 is an algebraic function of u; viz.: 

cos 6 = — B ± */p ± V q ± */r, (5) 

in which p, q and r are the roots of Euler's cubic, and the product of the radicals 

has the sign of — G, where G-= D — 3BC+ 2B S = — 4^- Since jb x and y x are 

negative, so are /? and 8 and consequently G. Hence either all three radicals 
in (5) must be positive or one positive and two negative ; also, none of the roots 
j>, q, r, can vanish. These latter are given by 

p=:B*-C+\V-J+*/M+\Z/— J— */W, 



q =&—0+fV-J+*/M+±V-J—*/M, (6) 

where in each case the radical denotes the principal root, and where 

1 /J 3 

w = — g H 7 r— , and — 27 M= P — 27 J 2 (the discriminant), 

JT and J being certain integral functions of the coefficients B, C, D and E. 

Now since cos $ is an algebraic function of u, it can be double-valued only 
by some radical having different signs on the two sides of the maximum normal. 
This requires that one of the outer radicals Vp, */q, or */r shall change sign ; 
for in the definitions (6) each sign has been definitely chosen. But, since G is 
negative always, this requires that two of the radicals shall change sign. More- 
over, in the ellipse cos 6 is a continuous function of the continuous variable u, 

*In this section the notation is that of Burnside and Panton in Theory of Equations, p. 131. 
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and as^>, q and r can not vanish, the two radicals which change must be always 
of opposite signs, and their difference must vanish for « = Wj. 

Now these two radicals are */q and Vr, and these are distinct from Vp 
for m = Mj, which follows from conditions for equal roots of Euler's cubic: 
p — 27 J" 2 = for two equal roots, and 7 = J = for three equal roots.* 
It has been seen that the first condition must be satisfied for u = u x ; (6) then 
shows that q and r become equal. That the latter conditions are not satisfied 
may be seen by a laborious direct computation of / or J, which may be avoided 
as follows : If p = q = r for u = u x , the quartic equation (4) will have a triple 
root, cos $ = cos \ as three of the four possible combinations of signs in (5) 
yield the same value. But while cos 6 = cos % is a double solution regardless 
of the position of the origin on the normal, the solution is triple only if the 
origin lie on the major axis or at the center of curvature for the given normal.f 
Both of these possibilities are excluded, since pj > p 2 > radius of curvature 
> normal length to a; -axis. Hence, for the positions of the center of force 
at present admitted, I'J^O for u = %. 

Finally, the proper signs being chosen, (5) becomes 

cos0 = — B — Vp± {*/q — */r), (7) 

where one sign is required when < < X, and the other sign when 2, < 6 < ^ . 

The next step of importance is to show that the difference V q — Vr 

contains the factor (u — Wj)*; i.e., Vg — «/r = V« — u t • 4* (u) where 

L 4* («) ^ 0, eo . For this purpose let \/q — */r be expressed as a fraction, 

thus : 

,— ,- V — 3 M 
A/q — Vr= 



(s/q + Vr) [(— /+ VM)i + 1 + (— J— VM)f] ' 

which shows, since every term in each factor of the denominator is positive at 
« = «!, that Vj — Vr contains the factor u — u x to the same power as does V M, 

* Bumslde and Panton, loc. cil., pp. 124, 144. 

t For the circle u = w, meets the ellipse in only two coincident points at («,, X), unless the circle be the 
osculating circle for that normal, in which case they meet in three coincident points. If the pole be upon OX, 
the double solutions cos 6 = cos A, and cosO = cos ( — X), make either a quadruple solution. 

Verification may be obtained by the test with the derived equations after expressing i, and y x as functions 
of X. (See a problem, by the author, in the American Mathematical Monthly, Jan. 1909.) 
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Since M is a polynomial (of degree 12) in u, vanishing for u=-u lt some 
integral power of u — 1% is a factor. To test directly the order of this factor 
is unnecessary, as it will appear in the next equation. 

Differentiate (7) with respect to u, denoting derivatives by primes; sub- 
stitute V -J+*/~M=R 1 , ^/ — J—^/~M—B i) ■^~ — ^T r = 2U > and 
1 j. 1 = 2_^ so that -i= ( L*= U-iV, and-^L — ^ = U+iV; 



Vq Vr -s/3' Vr \/ q ' s/r */q 

and re-arrange: 



(8) 



Now since the line = /I is a normal, -=^ = at (m,, X). The only denominator 
in (8) which vanishes at « = «! is a/J£ If if contains the factor {ti — u^", 
M' contains the factor (« — «i)" _1 ; hence, to have -Jt = for u = u lt it is 

necessary that — >i> — 1, or v<2. Therefore i>=l. Also, putting s/u—v^—x, 

inspection shows that x occurs to the power unity in U and in (^ — ^ J , and 

to the power zero in V, B', 2BB' — C, and (*+*\, (2BB'—C) may how- 
ever vanish identically). Whether J' carries a factor x is immaterial ; suppose 
x appears to the power (n — 2). 

Then equation (8) may be written in the form 

8in * 35 = \ {g + hx + hx% + lxH) > (9) 

where g, h, h and I are functions of u, which do not vanish at « = % unless 
identically zero. Moreover, from (7) one obtains 

sin 2 = 1 — [(JB + */pf =p 2 {B + Vp) (Vq- */r) + (Vq- VV) S ] 

— o -{■ CX — X , 
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where b and c are functions of u, not vanishing with u — u lt and not identically 
zero since = % :£ 0, — , when u = u x . Hence 

(du^ _ , b + cx — x* . x 

W0/ ~ to + Aaj + ^+Za;") 2 ' l ; 

The numerator of (11) contains an odd power of x. Even if the numerator is a 
factor of the denominator, the removal of this factor will still leave the same 
factor in the denominator, as the latter is a perfect square. And the conjugate 
factor needed for removing odd powers, if present at all, is squared. Therefore, 
(11) contains essentially an odd power of x. 

(du\ z 
-j-g) is a complex 

function of u; thus 

(W = * {U) = Fi (M) + * F * {U) > (1 2) 

where F x and F % are real for « ■< «i , but i F % (u) is real for u >■ u x . 



§ 5. Exclusion of Laws for Non- Axial Centers. 

The foregoing results permit a great restriction upon the location of the 
center of force within the ellipse. As is well known the force necessary for the 
description of any curve as a central orbit, the center of force being at the pole, 
is given at all points of the curve by 

d z u~* 



f= h '*' (» + §£)• 



where h is the constant of areas. There are an infinite number of laws according 
to any of which the force may vary to produce a given orbit with a given constant 
of areas and a given position of the center of force;* but by imposing hypo- 
thesis (I) at the start, a single law is obtained for the variation of the force along 

*Cf . F. L. Griffin : Certain Trajectories Common to Different Laws of Central Force, The Astronomical Journal, 
Vol. XXVI (1908), pp. 3-4. Inattention to this fact seems responsible for the use of hypotheses inadequate 
for the conclusion stated concerning binary systems by Professor Moulton [Introduction-to Celestial Mechanics, 
p. 80]. In this connection the writer desires to acknowledge his indebtedness to Professor Moulton for much 
helpful criticism of this paper, which had its origin in the detection of the error mentioned. 
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the orbit as a function of u alone. Thus, for the positions of the center of force 
considered above, the differentiation of (12) gives 

If the force varies according to this law throughout the plane, the resultant force 
is imaginary when u<^u lf that is, everywhere outside of the circle u = u x 
[or p = pi]- By the lemma above, real motion is impossible without that circle, 
so that hypothesis (III) is contradicted. This law and the corresponding positions 
of the center are, therefore, excluded ; that is, the center of force must lie upon 
a principal axis. 

[It may be observed in passing that the foregoing demonstration applies to 
exterior points and to those upon the ellipse as well as to interior points, and 

might have been used to exclude such points, if < a, •< — . Further it excludes 

points of the minor axis within the evolute (except the center), since from such 
a point there is a non-axial maximum normal. But in some ellipses the evolute 
does not contain the entire minor axis.] 

§ 6. Exclusion of Laws for Centers on the Minor Axis. 

Next suppose the center of force to be any point (0, d) of the minor axis 
Then, choosing that point as the pole and introducing polar coordinates, 
equation (2) becomes, since x 1 = 0, y 1 =id: 

e 2 p 2 sin 2 + 2dpsin0 + d 2 + (1 — e 2 ) (p 2 — a 2 ) = 0. (13) 



The solution of (13) involves V (1 — e z )(d z + a z e z — e 2 p 2 ), and so does the 

d z v 
dO z 



expression for the force (unless d= 0), as may be seen by forming -^p-. Now 



this radical is imaginary for p > - Vd 2 + a 2 e 2 , so that real motion is impossible 

outside a certain circle. Again there is a contradiction of hypothesis (III), 
restricting the center of force to lie upon the major axis. [It may be noted that 
points in the line of the minor axis exterior to, or upon, the ellipse are likewise 
excluded.] 
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§ 7. Laws Admitted by Hypotheses (I) and (III)' 

Finally let the center of force be any point (d, 0) of the major axis. 
[There is no loss of generality in restricting d to positive values, as this amounts 
merely to a choice of the positive direction of an axis.] Equation (2) becomes 
in this case 

e 2 p 2 cos 2 — 2dp(l — e 2 )cos0 + (1 — e 2 ) (a 2 — d 2 ) — p 2 = 0. (14) 

The solution of (14) is 

cos = — 2 Id (1 — e 2 ) ± Ve 2 p 2 + (d 2 — a 2 e 2 )(l— e 2 )]. (15) 

pe 

If d >ae, the radical does not vanish for real values of p. Its proper sign 
is ascertained by substituting the coordinates of a particular point, = 0, 

i = ■*('-'/> *<!-■"'> ce) 

e 2 (a — d) v ' 

which requires that the lower sign be chosen. Equation (15) may then be written 



ycos0 = att — V 8u 2 + y, (17) 

where 

a-d{\ — e 2 ), 8 = (d 2 — a 2 e 2 ){l-e 2 ), y = e 2 . (18) 



If d<iae, the radical vanishes; but the corresponding value of p, viz., p = -J — - 

is not taken in the ellipse unless d<ae 2 . Tests similar to (16) show that 

when ae 2 <d<Cae, the sign of the radical is always — , and 
when d<^ae 2 , the sign is + near = 0, and — near = n. 

[When d<^ae 2 there are two normals from the pole besides the axis; the sign 
changes at the foot of each.] 

Differentiate (17) twice with respect to 0, putting R — s/Bu 2 + y (positive): 

. / 8u\ d 2 u By /du\z _ 

-y coa6 = { a -^)dT 2 -R-^de) =~™ + R - 
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Substituting for (-tk) the value obtained by squaring and adding (17) and (19), 

one has 

/ Bu\ /d z u \ 8y iy z — (au — B) z _ 8u z y .. 

whence 

d*u y z (8 r + a*-8) 
W0 2 (aB — 8u) z ' K ' 

Multiplication by h z u z gives the required law of force [if d>ae 2 ] 

/ = ^ = 6 P (22) 

J (aR-8u) 3 (a^8 + yp z -8) s ' K } 



where 



e=h z y 3 (8 + - 1\ = W a z e 6 (1 - e 2 ) 2 . (23) 



[When d<^ae z , the corresponding law in any part of the ellipse is obtained 
from (22) by assigning to the radical a sign opposite to that specified above for 
the equation (17).] 

The laws represented by (22) have quite different properties according as 
d % ae; though in both cases / vanishes "at infinity". 

Type I, where d>ae. Here, since 8 >• 0, /is everywhere real ; also / 
is everywhere finite since the vanishing of the denominator would require 

p 8 = p| = v 8 -, a negative quantity; and/ is everywhere positive 

[attractive], for /3>aV , /^ + yp 2 would require p s < p| . 

Type II, where d<^ae. Here, since 8 < 0, / is imaginary within the circle 

Q 

o z = — -; also / becomes infinite at a finite distance p [only for the law 

r y 

required near = when d < a e % , — for the others as in Type I] ; and in the 
same case is negative [repulsive] for p>p . An exception is the very important 
case where a — : the center of force lies at the center of the ellipse and the 
law is that of the direct distance. 

Another special case of interest is that where d = ae z f the center of force 
being the center of curvature for the end of the major axis. The statement has 
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been made * that a particle can never, in a finite time, arrive at an apse, the 
center of whose osculating circle is at the center of force, but must approach the 
osculating circle asymptotically. The present example is a case to which the 
statement does not apply, since the time from any point (« , O ) to the apse 
(«!, 0) is given by 



«/«o (u 



where q> (%) dp. 0, oo . This integral is finite by a standard test. 

§ 8. Proof of the Operation of the Newton's Law. 

A property common to the laws of both types, except in special cases, is that 
the force vanishes at the center of force. For, from (22), it is evident that/= 
when p = 0, unless V (3 (± a — \//3) = which would require either d = ae 
or d = a. The latter exception has been excluded, since in that case the center 
of force would lie upon the ellipse. The former exception occurs when the center 
of force is at a focus; and in this case, since /3 = 0, (22) reduces to Newton's law, 
as it should. 

Therefore, as the vanishing of / for p = contradicts hypothesis (II), the 
Newtonian law must operate in the given binary system. 

[It is evident that the method of this section may be applied equally well 
to exterior points in the line of the major axis. And, again, the laws of type II, 
which prescribe an imaginary force in a part of the plane, might have been 
excluded otherwise by strengthening hypothesis (III) to require the possibility 
of real motion throughout the plane. Reasons for the choice of the method here 
adopted will be given in § 15.] 

Part II. Trajectories for the Laws of Type I. 

§ 9. Stability. 

Any law of types I and II admits the given ellipse as a central orbit, in 
which, at any instant, there are certain initial conditions. The question 
naturally arises as to the sort of curves that would be described if the initial 
conditions were different. 

* Cf. Routh, loc. eit., p. 283. This author probably intended, however, to have his statement understood 
as applying to only the special types of laws discussed in earlier pages, and not to be taken in the general form 
in which enunciated. 

10 
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With the exception of the law of the direct distance, whose trajectories are 
well known, no law of type II permits real motion in all parts of the plane. 
These laws will not be discussed further. All laws of type I (including that of 
Newton, whose simple trajectories will not be discussed here), produce for any 
real initial conditions real orbits, some of whose properties will now be 
established. 

Consider the law (22) for any chosen value of d, (ae<d<a). If a particle 
be projected from any point in any direction, with a velocity less than a certain 
function of the initial distance (the "velocity from infinity"), the orbit is stable. 
That is, the motion is periodic, there being both a pericenter and an apocenter ; 
and a slight change in the conditions of motion at any instant will produce only 
a slight change in the apsidal distances. 

At t=t let u = u , B = B ; and let the velocity be v and the constant 
of areas be h % . Then, if h z dp. 0, the differential equation of the path is, by (22), 

*K- + £9=p3=?ir (24 > 

Integration of (24) after multiplication by 2 -=^ gives 

* [■• + (ff )1 =*+m=*+ 'j%:%j$ <»> 

where Jc z is a constant, and where the left member equals the square of the 
velocity. It is easily seen that the limit of ^-f- F(u) as u becomes infinite is 

finite. Hence there is some value u x , such that, for u^>u 1} hi f -^ J < 0. 

But hi I -jTj- ) =f-jVJ. Therefore there is a circle u = u 1 , within which the 

radial velocity, — and likewise the motion, — is imaginary for the given initial 
conditions. This establishes the existence of the pericenter. 
From the initial conditions follows the equation 

h - vt - Ffa) = v$ -V \ ( 26) 

where V is the velocity acquired by a particle in falling from rest at an infinite 
distance to the initial point, and is real. The latter statement follows from the 
fact that e, a, y, (aB — /? w ), — and hence also 2aB — /? Mq, — are positive ; 
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so that F(u ) > 0. Now it is possible to choose v <^V ; and for any such choice 

&S5<0, and v 2 <0 at « = as follows from (25). Therefore h% (i£)* < ° 

for all values of u less than some value u 2 ; so that the motion is imaginary 
outside the circle u = Wg. Thus the orbit has an apocenter. 

Since in all central orbits included between two apsidal distances, the apsidal 
angle and time are constant, if the force be single-valued as required by (22), 
the radius vector is a periodic function of the longitude and of the time. 

If h % = 0, then -j- = and -^-s = 7 — w — ^— ^ . These two equations 
at at \ol H — p U) 

show the motion to be a periodic oscillation in the line = O . 

From (26) it follows that a small change in v would affect k z but slightly. 

And, since h z = v p sin 4- (where -^o is the angle between the radius vector and 

tangent at t = £ ), a small change in either v or 4 would produce only a small 

change in h z . Hence % and t%, the values of u for which -=-^ = 0, would 
change but little. The orbit is stable. 

§ 10. Other Oonios as Orbits. 

Another question presents itself: is the given ellipse a solitary conic orbit 
under this law of force, or are there others, — in particular, other ellipses ? 
It will be shown that there is a double infinitude of elliptic orbits, and that 
through every point of the plane there passes in every direction one, and only 
one, conic section which is a possible orbit. 

In the first place if any other orbit be an ellipse, the center of force must lie 
upon the major axis and not between the foci; for, if it were elsewhere, the law 
would be other than of type I. 

Hereafter let the constants of the original ellipse be denoted by the subscript 1, 
and the corresponding constants for any other elliptic orbit, by the same letters 
without subscripts. Then the equation of any such orbit is 



ycos(0 — v) = au— \f (3u z -f- y, (27) 

v being an arbitrary constant, the longitude of the pericenter. Also the law 
under consideration becomes in the new notation: 

f = (aaVfr+V^-fr) 8 ! (22l) 
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and similarly the force acting in the new ellipse is given by 

/= 6 P , (28) 

where e = h 2 a*e z (l — e s ) s , h being the new constant of areas. Hence for all 
values of o, the values of/ in (22 x ) and (28) are identical; from which it follows 
that the irrational terms in the two must be identical, and, likewise, the rational 

terms. The application of this principle to yj£- in (22 x ) and (28) gives 

« *SW _ «i V ft q */~y _ % V7i ^ _ ft c 2 q\ 

e * ^i— ' ~7W~- fil % ' e % — gj* * v ; 

From (29) it follows immediately that 

«_— «i_ §_ = §i L — ?k (30) 

Vy V/i' y 7i' y 3 7?' V ; 

That is, in all elliptic orbits which may be possible, the three quantities 

-2=, •£- and — ; must remain invariant. Let these be denoted by A, B and E, 
v y y y 8 

respectively; then a = Ae, fi = Be* and e = Ee s . The substitution of these 

values in (22 x ) gives for the law of force: 

f — E ? (31) 

r ~(AVB + 9 *—B) s ' K ' 

and in (27) gives 

ecos(0 — v) = Au— \/Bu? + 1; (32) 

so that any possible elliptic orbit under the law (31) is represented by (32) for 
some value of e, — the original ellipse for e = e x . 

Conversely, it will be shown that for all values of e and v the curves (32) 
are conies, — a double infinitude possessing certain properties because of the 
invariants A and B. And further there are real initial conditions, — determined 
by two invariants, E and another (obtained later), — for which each member of 
the family of conies (32) is a possible orbit under the law (31). 

§ 11. The Family of Conies. 

Consider the equation (32) for arbitrary values of e and v. All the distinct 
curves which are represented by it are obtained for 0<e< oo and O^v^ln, 



Griffin: On the Law of Gravitation in the Binary Systerns. 77 

the curves obtained for negative values of e being included in these. For, by 
changing the sign of e and adding ± n to v, the equation is unchanged. 

To vary v while keeping e constant does not change the shape of a curve, 
but merely rotates it in the plane. It is, therefore, sufficient to consider the 
curves for which v = 0. This particular infinitude will hereafter be called the 
restricted family ; its geometrical properties are possessed by all families in which 
v has a fixed value. 

If rectangular coordinates be introduced by cos — x u, sin = y u, the 
equation (32) with v = becomes: 

x z (l — e z ) + 2Aex + y z = D, (33) 

where 

D = A z — B = (af — d\ ) (1 — e\) > 0. 

Equation (33) represents for any value of e a conic section, whose transverse 
axis is in the aj-axis, and whose eccentricity is e. For e dp. 1, the center is 

( — 2 , 0j, and the square of the major semi-axis is a z = . e "*" 2 2 . The 

apocentral distance in any ellipse is, then, 

(Be z + D)* + Ae D 

1-e* ° r (Be z + D)i-Ae ] 

similarly the pericentral distance in any conic is , R 2 4 . . The former 

steadily increases with e, having its smallest value VD for e = 0, and becomes 

infinite as e == 1. The latter is a decreasing function of e, having its largest 

value VJ5 for e = 0, and approaches zero as <?==<» . Hence every conic of the 

general family has its pericenter within, and its apocenter without, the circle of 

the family, p z = D. 

Various geometrical properties, belonging chiefly to the restricted family, 

are readily seen. As neither of the invariants A and B is so well suited to inter- 

A 
pretation as two of their combinations, viz., D and — ^=, these latter two will 

be used. 

I. The First Invariant, D. All conies of the family with collmear axes pass 
through two fixed points; viz., the points where the circle o z — D meets the 
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perpendicular to the axial line through the origin. For equation (32) with v 

fixed is satisfied by p = «/Z>, = v ± - , whatever be the value of e. [Another 

proof rests upon the fact that all members of the restricted family pass through 
(0, ± VZ>).] 

It is possible to define new families of conies related to the family (32) in 
which there is a simple interpretation of this and the following invariants. If a 
new conic be described, having the same eccentricity as some given conic of the 
restricted family, and having as its transverse axis the segment between the 
center of the given conic and the pole with respect to the given conic of the line 
joining the fixed points (0, ± V2>), then this new conic passes through the same 
two fixed points, whatever the value of e* To establish the property it is sufficient 
to note the coordinates of the vertices of the new conic, 



(- r 4^,0)and(^,0) 



Evidently the auxiliary circle of the new conic cuts the F-axis where F 2 = 



1-e 2 ' 

F being the mean proportional between the two segments of the axis. Let y be 
the ordinate in the new conic at a; = 0. Now, in any ellipse, y z /Y z = b 2 /a 2 ; 
and, in any hyperpola, y 2 / F 2 = — b z fa z . In either case, y z = F 2 (1 — e 2 ) = D. 

II. The Second Invariant, — == . In all conies of the restricted family, the 

slope at either fixed point bears to the slope at the end of the latus rectum a ratio 

which is constant. For the latter slope is ± e, while the former is seen from (33) 

Ae A 

to be ± — r -. The ratio is ± — ,==, which is constant, 
v D v D 

If a new conic be described, having the same center and eccentricity as 

some given conic of the restricted family, and having a focus at the origin, the 

new conic passes through two fixed points whatever the value of e, so that its 

latus rectum bears a constant ratio to the segment between the fixed points of the given 

family. For, the major semi-axis of the new conic being — , its latus rectum is 

e 

* The members of one family do not belong to the other (except for a single value of «). For in each 
family there is but one conic with a given eccentricity, and corresponding conies have axes of different lengths. 
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2 d (1 e 2 ) 

or 2 A ; and the distance between the fixed points (0, =L V D) being 

^ 

2 s/ D f the ratio is .-=, • Again, since the constant latus rectum lies always in 

the F-axis, all the new conies pass through (0, ± A). 

III. An Alternative Invariant, B. The family defined in II intersect the 

tangents at the vertices of corresponding members of the given family on two straight 

lines, y 2 = B; and conversely, on y 2 = — B. These statements may be proved 

analytically by using the equations of the tangents, or geometrically by observing 

that the segments into which the major axis of the new conic is divided by the 

d d 

vertex of the old are a and - + a. It follows that the tangent is cut by the 

e e 

d® /r^ g^ 

auxiliary circle of the new conic where y 2 = 2 , and hence by the conic 

(d% a 2 e 2 )(l e 2 ) 

where y 2 = * /-^ = B. Similarly for the converse. 

IV. Any two conies of the general family intersect, so that no one lies wholly 
within any other. For let the values of v and e in the two curves be respectively 
v 1} e x and v s , e z ; and let the radii vectores be denoted by o x and p a . Suppose 
e 2 > e 1 . Then, by the remark above concerning the pericentral distance, 
a % — d % < a x — d\ ; and hence, when = v % , p 2 •< pi , f° r a i — <h < Pi • 

Case I. When e 8 < 1, both orbits have apocenters, and a z -\- d 2 > a t + d\ . 
Hence, for 6 = v z + n, p 2 > p x ; and, for some value of 6 between v z and v 2 + n, 
the two ellipses must intersect. 

Case II. When e 2 > 1 >■ e 1 , the ellipse has an apocenter, while the other 
conic extends to infinity. Hence there are values of 6 where p 2 >pi, so that 
the curves cut. 

Case III. When 6i>l, there exist values of 6 where p 2 >pi. For the 
value of 6 — v where any conic of (32) goes to infinity is <£> = arc cos f J , 

showing that $ > - and $ is a decreasing function of e. Now p is finite from 

d = v — $ to d=zv-{-q>; and since $i !> # 2 > ^c , the range of values of for 

which p 2 is finite is smaller than that for which p x is finite, but is greater than n. 
Hence at least one of the radii vectores where p 2 becomes infinite lies in the 
region where p x is finite. Near this, p 2 >pi. 
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V. Through every point there passes in every direction one, and only one, conic 
of the family (32). For, let the point be,(u , O ), and let the given direction 
make with the radius vector to that point the angle i/> . Then, in the required 

/<*0' 



curve, tan ^o = — w o(t— ) > where from (32) : 



-§A— A*/ B*u 2 +\—Bu r34) 

du VBu 2 +l 2 — (Au — V^Bw 2 + l) a ] v *■ ' 

Hence if any member of the family is to have 4* = 4*0 at m = « , its eccentricity 
must have the positive (or zero) value given by 

e 2 = (Au - V B ul+iy + ^A-^^^ful cot 2 ^ . (35) 

Evidently el > 0, so that there is a real positive (or zero) solution. The sub- 
stitution of e = e z , 6 = d and u = u in (32) determines v so that the conic will 
pass through the given point. The equations determine both c 3 and v uniquely, 
subject to the original restriction. 

Corollary. Every point in the plane is an apse in some conic of the family. 

For then 4o = k , and equations (35) and (32) become e z = ± (Au — */Bul + 1), 

and cos (0 O — v) = db 1, where the upper or lower sign is to be taken according 
as A*>B + ol, — or ? 1<D. Further, if pg<Z>, d =v; while, if o%>D, 
d = v ± n. That is, any given point within the circle f = D is the pericenter 
in some conic of the family ; and any point without that circle is the apocenter 
in some conic of the family. [This corrollary establishes the converse of a 
property obtained above.] 

§ 12. The Initial Conditions. 

It has been shown that, if the law (31) admits as trajectories any ellipses 
other than the original orbit, those ellipses must be members of the family of 
conies just discussed; and it is further necessary that the constants of areas have 
values determined as follows : 

The Third Invariant, E, or h z a 2 (I — e z ) z . The constant of areas must be 
proportional to the curvature at an apse. For the radius of curvature at an apse 
is a (1 — e 2 ), so that the constant of areas h divided by the curvature is invariant. 
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Several questions now require consideration. 

I. Is an arbitrarily chosen constant of areas h % admissible for some conic of the 
family ? If so, then between h % and the required eccentricity e 2 must exist the 

relation E~K\ (& + a *~ @n) = h% {Bel + D), which requires 

h\Be\ = E—h\D. (36) 

Equation (36) defines a unique positive number e a , provided simply that ^<^ ( 

the value of h* in the circle of the family. [If h % = 0, e s = <x> , the curve being 
a straight line, as is indeed evident from other considerations.] 

This result does not show, however, that there is an upper limit to the initial 
velocity of projection. For h = vo sin 4* constantly; and therefore arbitrarily 
large values of v yield admissible values of h, if only the initial distance or the 
angle 4* be sufficiently small. 

II. For an arbitrarily chosen constant of areas h % , less than the limit mentioned 
above, what further restriction must be put upon the initial conditions in order that 
the trajectory be a conic? Let the values of o, u, 0, v and 4- at the initial point 
be designated by the subscript zero. Now in order that the trajectory be a conic, 
e must have the value e a , defined by (36). 

Hence, in the first place, p must lie between a % — d % and a 2 + d 2 [or oo 
if e z > l]. Further 4o must equal the value of 4 at the distance p in the conic 
for which e = e 2 . And, finally, v must be given by w = h% u esc 4'o- 

To formulate the second of these conditions, consider equation (26). Since, 

dd /du\ z 

in any curve, tan 4 = — u j~> or ^csc 2 ^ = « 2 + ( -rs J , the value of wgesc^o 

in the orbit (25) is, for arbitrary conditions of projection, 



u esc 4, - M + A z hKAVBul+1 _ Buor (37) 



And in the conic for e = e s , the value at u = u, 



is 



u* esc 2 J, - u* 4- (^^ + D[^-(^»o-^^^+l) 8 ] , 38 ^ 

Equating the right members of (37) and (38), and observing that 

E = hHB{e»-l) + A*l 
11 
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h 
a condition on j4 is obtained. The irrational parts must be identical, likewise 

the rational parts: 

Ka 



£f(— 2A 8 Bu )V Bu 2 + 1 = — 2ABu (ei — l)V Bui + 1 

ft < 39) 

U A* [A* + « S (A* B + B*)-] = (el — 1) [A 2 + u§ (4 a 5 + J5»)]. 

"2 

Each equation requires that vf = g . a ; that is, for all orbits of the family, 

— ±-r is a fourth invariant, restricting the initial conditions. 

For a conic orbit with the constant of areas h z , the velocity of projection 
at the distance p must be » , where 

v% = ft 2 + F(u ) = hi e -t^± + F(u' ); (40) 

and 1^0 is given by h % = v p sin 4v Or, reversing the order, ^o ma y be 
obtained from (38), and then v from either of the two equations last written. 

III. To determine the velocity, with which a particle projected at any point in 
any direction will describe a conic. There is one, and only one, conic through the 
given point in the given direction. Since « and ^o are known, (35) determines 
e % , real and positive; and (36) then determines h % . Finally v is given by (40) 
or by v = h 2 u esc i^ . It is evident that there is one, and only one, real positive 
value of » for which the orbit is a conic. The fact that there are initial velocities 
for which the orbits are other curves than conies, is important in a later section. 

§ 13. The Apsidal Angle in the Non-Oonic Orbits. 

The question as to whether "the line of apsides advances or recedes," 
— that is, as to whether the apsidal angle © is greater than or less than n, — 
is a matter of interest. The expression for the force in laws of type I may be 
written, / = u z P («), where 

Now it has been shown elsewhere* that, if P' (w) > for u t > u > « 2 , then 

* F. L. Griffin : On the Apsidal Angle in Central Orbits, Bulletin, of the American Mathematical Society, 
Vol. XIV (1907), pp. 6-16. 
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© >?t in all orbits lying wholly between the circles u — u x and u = Wg; and con- 
versely. Similarly, if P'(«)<0, © < n. And further, if P' («)< 0, and 
decreases with u or is constant, then in any family of orbits, one of whose 
apsidal distances is constant, © increases with the second apsidal distance. 
Analogous statements apply in the opposite case. 
For type I, 

Hence P'(«)>0 if B>Au; that is, if p a >Z>; while P'(w)<0 if p 2 <l>. 
Consequently in all the orbits which lie wholly within the circle p 2 = D, © < n ; 
and in all orbits wholly without the circle p 2 = D, © > n. [It is interesting to 
compare this result with the fact that in the conic orbits © = n, and to note that 
a new proof is thus obtained for the earlier statement that every conic orbit 
admitted by the law (31) has its pericenter within and its apocenter without the 
circle p 2 = Z>.] 

Moreover, since P' ( J = and P' («)<() when u «< — , there 

exists a region within the circle p a = D and without another concentric circle, 
throughout which © decreases as the apsidal distances in the orbits decrease 
The corresponding question for orbits outside p 2 = D requires more extended 
computation. 

§ 14. A New Proof with Classic Hypotheses. 

From the foregoing discussion, it is possible to obtain a new proof of the 
operation of the Newtonian law in the binary systems, based upon the hypotheses 
mentioned in § 1 in connection with the classic theorem of Halphen and Darboux. 

The classic hypotheses imply those of the present paper ; for (I) and (II) 
are common, and, if every orbit is a conic, then the law according to which the 
force varies (not merely in one orbit but throughout the plane) is such as to 
permit real motion everywhere. And, since the operation of Newton's law 
follows from the new hypotheses, it follows also from the old ones. 

Again, a new proof is possible with the assumptions of Halphen and Darboux 
modified as suggested in § 1. For the second theorem of Bertrand referred to 
may be applied, with assumption (III) and the fact that at least one conic orbit 
is an ellipse, to show that the force is central. Consider that ellipse ; the center 
of force may lie only upon the major axis. For, if it were elsewhere, the use 
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of (I) would furnish a law, for which some initial conditions would give rise to 
imaginary orbits, [§ 6], in violation of (III). Moreover, if the center of force 
were any point of the major axis save the center [excluded by assumption] or a 
focus, the corresponding law would admit imaginary orbits (if of type II) or real 
non-conic orbits (if of type I), [§ 12], again in violation of (III). Thus, New- 
ton's law must operate. 

§15. Possible Modifications in Hypotheses. 

(1) As observed above, laws of type II might be excluded by strengthening 
hypothesis (III) to require the possibility of real motion, not merely sufficiently 
far from the center of force, but throughout the plane, an increase in severity 
which may seem on first consideration to be extremely small. 

It is true, however, that the hypothesis in this revised form would exclude 
so important a law as that according to which an oblate spheroid attracts an 
exterior particle in its equatorial plane ; for, if the latter law be required to 
operate within the spheroid, the force is imaginary within a certain interior 
circle.* But it is well to observe that even this exclusion would not be excessive 
in the present problem, inasmuch as the latter law admits as orbits no conies 
except circles. 

Nevertheless, if the attracting body be regarded as having size, the severity 
of the revised hypothesis is greatly enhanced, as it would exclude the operation 
at exterior points of any law which, if operating also at interior points, would 
prescribe imaginary values of the force within the attracting body, even though 
all (physically possible) real initial conditions give rise to real orbits. If 
modified to require only the reality of the force everywhere outside the attract- 
ing body, the hypothesis would not suffice to exclude laws of type II for which 
ae > d > ae 2 , and would suffice to exclude those where d < ae 2 only when 
strengthened by the assumption that the bedies do not collide. [These latter 
statements rest upon the fact that, if pi denote the value of p for which the 
imaginary radical vanishes, p a < a — d (the least radius vector in the orbit) if 
ae^> d^> ae 2 , and p 2 = the minimum normal if d < ae 2 .] When ae^>d^> ae 2 , 
however, the center of force would at least have to be sufficiently near the focus 
to have the principal star contain the focus ; for the inequality p 2 > ae — d 
follows from the known inequality > (d — ae) (d + ae — 2ae 8 ). 

*For the law in question see Moulton. loe. cit., p. 133, or paper cited in Bulletin, p. 13. 



Griffin : On the Law of Gravitation in the Binary Systems. 85 

(2) Again, it has been noted that all exterior points, and all points on the 
ellipse except the ends of the major axis, could be excluded without the aid of 
the reference cited in § 1. These two, and all other, points of the curve might 
be excluded by formulating as an assumption for material points the physical 
principle of impenetrability, or by otherwise prohibiting collision. This change 
would, however, involve an altogether unnecessary (even if mild) addition to 
the assumptions of the paper. 

(3) Finally, in any non-rigid body tides would be formed with a resulting 
difference in the attracting force, which, considered as a perturbation, would in 

general forbid the constancy of v , (the fourth invariant). Doubtless a 

substitute for hypothesis (II) can be framed, based upon some such requirement 
as that, under tidal disturbance the force shall cause the velocity at any instant 
to be such as to permit the description of a conic as the subsequent orbit were 
tidal disturbance to disappear. In addition, it would be necessary to require 
that the center of force be not at the center of the orbit, in order to exclude the 
law of the direct distance which compels an elliptic orbit for all initial conditions. 

Williams College, September, 1908. 



